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In the last years many studies have been devoted to the transition from 
classical mechanics to quantum mechanics. These studies are motivated by 
the so-called quantum chaos (see Ozorio de Almeida 1990, Gutzwiller 1990, 
Casati and Chirikov 1995). An important aspect is the semiclassical quan- 
tization formula of the energy levels for integrable and quasi-integrable sys- 
tems, i.e. the torus quantization initiated by Einstein (1917) and completed 
by Maslov (1972, 1981). As is well known, the torus quantization is just 
the first term of a certain ^-expansion, the so-called WKB expansion, whose 
higher terms can be calculated with a recursion formula at least for one degree 
systems (Dunham 1932, Bender, Olaussen and Wang 1977, Voros 1983). 

Recently it has been observed by Prosen and Robnik (1993) and also 
Graffi, Manfredi and Salasnich (1994) that the leading-order semiclassical 
approximation fails to predict the individual energy levels within a vanishing 
fraction of the mean energy level spacing. This result has been shown to 
be true also for the leading (torus) semiclassical approximation by Salasnich 
and Robnik (1996). 

In this paper we analyze a simple one-dimensional system for which we 
are able to perform a systematic WKB expansion to all orders resulting in 
a convergent series whose sum is identical to the exact spectrum. For this 
system we show that any finite order WKB (semiclassical) approximation 
fails to predict the individual energy levels within a vanishing fraction of the 
mean energy level spacing. 

The Hamiltonian of the system is given by 

H=£- + V(x), (1) 

where 

Vtx) = ■ (2) 

cos (ax) 

Of course, the Hamiltonian is a constant of motion, whose value is equal 
to the total energy E. To perform the torus quantization it is necessary to 
introduce the action variable 

1 r \/9m i — 

/ = I pdx = ^(VE - Ju ) . (3) 

The Hamiltonian as a function of the action reads 
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and after the torus quantization 



I = (v + \)h, (5) 



where v — 0, 1, 2, . . ., the energy spectrum is given by 



EF = A[( V + \) + \b\\ (6) 

where A = a 2 h 2 /(2m) and B = ^8mUo/(ah). 
The Schrodinger equation of the system 

l-^^ 2 +V(x)Mx) = E^x), (7) 

can be solved analytically (as shown in Landau and Lifshitz 1973, Fliigge 
1971) and the exact energy spectrum is: 

El* = A[{v+\)+ l -Vl^W]\ (8) 

where v = 0, 1, 2, . . .. We see that the torus quantization does not give the 
correct energy spectrum, but it is well known that the torus quantization is 
just the first term of the WKB expansion. To calculate all the terms of the 
WKB expansion we observe that the wave function can always be written as 

i 

ijj(x) = exp (-o-(rr)) , (9) 

where the phase a(x) is a complex function that satisfies the differential 
equation 

a'\x) + (-)a"(x) = 2(E - V(x)) . (10) 
% 

The WKB expansion for the phase is given by 

fc=0 1 

Substituting (11) into (10) and comparing like powers of H gives the recursion 
relation (n > 0) 

a'l = 2(E - V(x)), jr a' k a' n _ k + a" = . (12) 

fc=0 



3 



The quantization condition is obtained by requiring the uniqueness of the 
wave function 

/°° h r 
da = Y,(~) k f d(y k = 2nhv , (13) 

where v — 0, 1, 2, . . . is the quantum number. 

The zero order term, which gives the Bohr-Sommerfeld formula, is given 

by 

p p / TP 1 

jda = 2j dx\j2(E - V{x)) = 2vr^(W - - -B) , (14) 

and the first odd term in the series gives the Maslov corrections (Maslov 
index is equal to 2) 

(-) j> d<7i = (-)- lnp\ contour = -nh . (15) 

The zero and first order terms give the equation (6), which is the torus 
quantization formula for the energy levels (Bohr-Sommerfeld-Maslov). Here 
we want to analyze the quantum corrections to this formula. We observe 
that all the other odd terms vanish when integrated along the closed contour 
because they are exact differentials (Bender, Olaussen and Wang 1977). So 
the quantization condition (13) can be written 

00 fi r 1 

E(-) 2 "f da 2k = 2irh(u+-), (16) 

thus again a sum over even-numbered terms only. The next two non-zero 
terms are (Narimanov 1995, Bender, Olaussen and Wang 1977, Robnik and 
Salasnich 1996) 

(-Y & da 2 = = — —— / dx . v = , (17) 

V J y/2^l2dE*J ^ E _ V{x) 1 } 

,h, A / , ft 4 . i d 3 r , v" 2 (x) i d 4 r , v 2 {x)v"(x), 

(—) f da 4 = ^rj-[ 7— - / dx 7—— / dx- 



, j (2m)3/2 L 120 dE*J ^J E _ V{x) 288 dE± J ^E-V(x)' 

(18) 

A straightforward calculation of these terms gives (see the Appendix) 

( )./ (fal= (19) 



4 



and 

r h sA r , 2nh 



Up to the fourth order in H ~ £> the quantization condition reads 

*f = ^ + 5) + ^ + ^-^. PI) 

The first two terms on the right side give the torus quantization formula, 
and the other two terms are quantum corrections. Higher-order quantum 
corrections quickly increase in complexity but in this specific case they can 
be calculated. We first verify by induction, following Bender, Olaussen and 
Wang (1977), that the solution to (12) has the general form 

a'jx) = (^) 1 - 3n P n (cos (ax)) sin /(n) (ax) , (22) 

where f(n) — for n even and f(n) = 1 for n odd, and P n is a polynomial 
given by 

s(«) 

P„(cos (ax)) = C n ,i cos 2 *- 3 " (ax) , (23) 
1=0 

with g(n) = (3n — 2)/2 for n even and g(n) = (3n — 3)/2 for n odd. 

The integrals in (16) arc performed by substituting z = tan (ax). In this 

way the 2/c-term reduces to 



(24) 

We observe that 

(1 + Z 2)3fc-I-1 



(E-U - Uoz 2 )™- 1 / 2 



^-1 r (|) a 3 *" 1 r, d + z 2 ) 8 *- 1 - 1 



( 1} " r(3fc- 1) di^- 1 f dz (E-u - u z 2 y/ 2 ' (25) 

so the only non-zero term is for I = 



d 3k-i ( 1 + z 2yk-i 2 6k ' 1 T(3k- 1/2) 2 9 



,3k- 1 



//r 1 I t — 1 : : R 



3fc-l 



dE™- 1 J (E — U — Uqz 2 ) 1 / 2 u 1 ! 2 r(6fc-l) <9£ 3fc 
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= r( 6 fe - 1) r ^^f^ 27r ' ^ 

where (5 — (E — U )/U . At this stage we obtain 

(^ k fda 2k = . (27) 

Now we need to find the coefficient C 2 k,o explicitly. By inserting (22) with 
(23) in the recursion relation (12) we obtain 

n n—1 

J2 C k,oC n „ kfi -(2mU a)C n ^ = J2 C kfi C n „ kfi +2C nfi -(2mU a)C n ^ lfi = , 

k=0 k=l 

(28) 

from which we have 

Y k-l 

C kfi = -[(2maU )C k - lfi - £ C jja C k -j, Q ], C , = 1 . (29) 

From this equation one shows Ci )0 = maU . Further, it easy to show that all 

higher odd coefficients vanish, C2fc+3,o — for k — 0, 1, 2, The solution of 

this equation for the remaining nonzero even coefficients is given by 

C 2kfl =(-l) k (2mU a) 2k 2- 2k (^j , (30) 

which can be verified by direct substitution in equation (29) resulting in an 
identity for half integer binomial coefficients. Then the integral (27) can be 
written 

(5)»/,fo* = (-l)fi-2 TO -'(2 m )" 2 -'2-(|)t/ '-" 2 = . 

(31) 

In conclusion, the WKB quantization to all orders (16) is 

E ^- A ^ + \> + \W)^ f - (32) 

Because Y,T=q (j) B l ~ 2k = y/1 + B 2 we have E™ = E^°\ i.e. the WKB 
series converges to the exact result (8). 
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Now we can calculate the error in units of the mean level spacing AE U = 
— Ef, x between the exact level Ef, x and its WKB approximation E^ 
to Nth order: 



E ex — E^ 



AE. 



V 




for v — > oo . (33) 



The limit clearly shows that even for arbitrarily small but finite ft (1 < < B < 
oo), the relative error for any finite WKB approximation becomes constant, 
and scales as 



rpex pp(N) 

'^^lv 2 Jvr^rr B^oc. (31) 




Note that the limit B — > oo is equivalent to the limit ft — > 0. For our present 
system we can conclude that to any finite order semiclassical approximation 
the error measured in units of the mean level spacing remains constant even if 
the quantum number increases indefinitely, contrary to the naive expectation. 
This confirms the general statements made by Prosen and Robnik (1993). 

Our approach leading to the above conclusion rests upon a systematic 
WKB expansion for the potential V(x) — Uq/ cos 2 (ax) using the technique 
of Bender, Olaussen and Wang (1977). We are able to calculate all orders, 
the series is convergent and can be summed precisely to the exact result. 
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Appendix 



In this appendix we show how to obtain the formulas (19) and (20). In all 
integrals of this section the limits of integration are between the two turning 
points. After substitution z = tan (ax), we have 

V'\x) AaU 2 ryfi z\z 2 + 1) 



f dx v _ m = r 



dz- 



where (3 — (E — U )/U . In conclusion we have 



Ui %f + 4P)1), (35) 



: i»l fc = -s« = -¥' (36) 



with B = y/8mU /(ah). 

To obtain the formula (21) we proceed in the same way. 

fdx V " 2 ^ 4a 3 £/ 2 ryfi ^ {9z 4 + 6z 2 + l){z 2 + 1) 

' x ^e- v(x) VUo~ z VP-? 2 

= -^(45/3 3 + 90/3 2 + 56/3 + 16)^, (37) 

From which we obtain 

m J iix ^wr ■■ (38) 

For the last integral we have 

V' 2 (x)V"(x) 8a 3 U 2 ry/P J z 2 (3z 2 + l)(z 2 + I) 2 

dz- 



J Jf, - VOA ^/C^ J- 



z 2 



SaU jl(105p 4 + 280/3 3 + 240/3 2 + 60/3) , (39) 



'U y ' 128 
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from which we obtain 
d 4 



In conclusion we have 



<7>' 



doA 



V ,2 {x)V"{x) _ 3157m 2 



dx 



dE'J p - V(x) 



h 4 



a 



TT a , 1 135 



J J 2 



(2m) 3 / 2 
64(2m) 3 /2[/2 



120 2 

2nh 



165 3 
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